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We analyze the scaling laws for a set of two different species of long flexible polymer chains 
joined together at one of their extremities (copolymer stars) in space dimension D = 2. We use a 
formerly constructed field-theoretic description and compare our perturbative results for the scaling 
exponents with recent conjectures for exact conformal scaling dimensions derived by a conformal 
invariance technique in the context of D = 2 quantum gravity. A simple MC simulation brings about 
reasonable agreement with both approaches. We analyse the remarkable multifractal properties of 
the spectrum of scaling exponents. 



Polymers serve as a testing ground for a field theory 
since the early seventies |l], ^ . The path integral formu- 
lation H allows for a direct interpretation of the paths 
as conformations of random walks (RWs) and self avoid- 
ing walks (SAWs). Closed paths in this theory have a 
multiplicity of the number of components M of the field. 
The formal limit M = excludes these loops and yields 
the polymer limit of field theory, that describes self and 
mutually interacting paths. The field theoretical descrip- 
tion of the scaling properties of polymer chains in good 
solvents has been generalized to the case of multicom- 
ponent polymer solutions Q and linked polymers [|. 
Recently, the theory was suited to describe polymer net- 
works of different species |fl 0. The simplest nontrivial 
case of a heterogeneous polymer network - a star-shaped 
copolymer in D = 2 dimensions - is the subject of the 
present report. 

Polymer field theory usually has to be evaluated in 
terms of a truncated perturbation theory series. Obvi- 
ously, it is always of special importance to compare these 
with the exact results if such exist. As a rule they are 
available for D = 1,2. For the critical exponents that 
describe the scaling of homogeneous polymer chains and 
star polymers (or, more generally, homogeneous polymer 
networks) the perturbative approach S is in fair agree- 
ment with the exact data for D — 2 . For the heteroge- 
neous case however, only perturbative results |t], |j| were 
available until recently. An exact solution for scaling ex- 
ponents has been recently proposed that uses methods 
of conformal invariance on random graphs j^]. Below 
we undertake a comparative analysis of the properties 
of the perturbative solution and of the exact one and 
complement these results by MC simulations. We show 
where the approaches agree and explain why they differ 
in general. In this respect our paper confirms some of 



the speculative results of || derived from the algebra of 
conformal dimensions for copolymer stars on 2D random 
graphs (quantum gravity). 

A remarkable feature of the scaling spectrum of a 
copolymer star is that it possesses multifractal properties 
that can also be analyzed in terms of field theory and the 
conformal invariance approach. |?], ||, ||, [n], [l^] . 

a. Copolymer stars and their scaling exponents Star 
polymers as the most simple nontrivial examples of poly- 
mer networks || ^| are created by linking together the 
end points of polymer chains at a common core. Let us 
now consider a general star polymer made of chains of 
two different species (while all chains have equilibrium 
size R). When such a copolymer star is immersed in a 
good solvent its asymptotic properties are universal in 
the limit of long chains |0, j|. In particular, the num- 
ber of configurations (the partition function Z*) scales 
with R on some scale I as: 
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where A/, /, constitute a family of copolymer star ex- 
ponents n, Kj. These are universal and depend only on 
the space dimension D and the number of chains of dif- 
ferent species 1 and 2 (/i and fa correspondingly). If 
a nonvanishing mutual avoidance interaction is present 
between the chains of species 1 and 2 polymer field the- 
ory predicts two nontrivial physically different regimes 
Q for the scaling behavior of true copolymer stars de- 
scribed by Eq. (|l|): In the first case the polymers of 
species 1 behave as RWs and the polymers of species 2 
as SAWs, whereas in the second case both species display 
RW behavior jl3]]. For these two relevant cases a renor- 
malization group description of the exponents Xf t j 2 nas 
been given in terms of the scaling dimensions of appro- 
priately composed composite operators of polymer field 
theory In particular, the resulting perturbation the- 
ory was derived in the form of an e = 4 — D expansion 
as well as evaluated directly for fixed D Q. Follow- 
ing Ref. [Q where the copolymer star exponents rjf l j 2 
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were calculated, the e-expansion for the A-exponents of a 
copolymer star composed of two mutually avoiding sets of 
RWs (A? f 2 = —Vfrji) an d f° r a star of mutually avoid- 
ing sets of SAWs and RWs (A^ = -r?/ lj2 + fm¥,o) 
read: 



= - (/a - 3 + A) J[l - (/i + 

/ 2 +3C(3)-3)|]}, 

A / 1 ,/ 2 ( £ ) = -/i( 3 - h ~ 3/ 2 )| - /i(43 - 33/! + 

8/f - 91/ 2 + 42/i/a + 18/|)^ - /i[675 - 969/i 

256 

456/i 2 - 64A 3 - 2463/ 2 + 2290/i/ 2 - 492/^/2 + 
1050/ 2 2 - 504/!/ 2 2 - 108/ 2 3 - C(3)(712 - 936/i 4 

224/x 2 - 2652/ 2 + 1188/i/ 2 + 540/ 2 2 )] 
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where £(3) ~ 1.202 is the value of the Riemann zeta func- 
tion. While we do not display the exponents calculated 
at fixed space dimension D here explicitly, we show their 
numerical values below. 

The pcrturbative formulas of the kind given in Eqs. 
@-(||) may be evaluated for D — 3 where the topo- 
logical complexity of the problem does not allow for an 
exact treatment. One should be aware that for similar 
reasons an exact expression for the partition function Z* , 
formula (Q), can not be obtained for D = 2 either. How- 
ever, appropriate identification of the universality classes 
of D = 2 critical behavior with conformally invariant the- 
ories often allows to extract the exact values of critical 
exponents. In Refs. |j) this was successfully performed 
for the D = 2 copolymer star of two or more mutually 
avoiding bunches of SAWs and RWs. Here, we analyze a 
copolymer star of one or two RWs or SAWs of species 1 
and a bunch of fa RWs of species 2. In the notation of 
the above equations the following exact expressions for 
the corresponding exponents are derived from Refs. M: 
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with a = V24/ 2 + 1. 

As far as this is not an exact solution in terms of the 
original theory but rather an "exact conjecture" for the 
exponents we performed one more check of the scaling ex- 
ponents by a MC simulation for the simplest case when 
the copolymer star consists of RWs only (^ff 2 , ^2/2 )• ^ n 
these simulations the stars are grown on a square lattice 
until a non- allowed intersection occurs. The number of 
stars C (N) generated during the growing process is accu- 
mulated for all chain lengths N. Exponents are extracted 
assuming a power law C(N) ~ j\f~"° x fi- f 2 where the RW 
correlation length exponent is vq = 1/2. We generate 
stars up to a maximum chain length N maK = 10 3 and 
accumulate configurations until C(A max ) = 10 5 . The 



number of successful attempts to grow stars with longer 
chains decreases rapidly for higher fa + fa which increases 
simulation time drastically. While for the case 1 + 1 a 
total of 10 6 attempts were needed, this number rose to 
5 • 10 9 for 5 + 1. Therefore, we report results only for 
fa + fa < 6. From the bare simulation data we extracted 
the numerical values of the exponents by an extrapola- 
tion to l/N — > 0. In Table | we report the results and the 
statistical error of this extrapolation in the lines marked 
'MC and '±A\ We compare with the exact numbers 
according to Eqs. (Q). The results are in fair agreement 
with a slightly growing discrepancy as the number of RWs 
fa increases. 

In Fig. [l] we compare the results with the perturbative 
data of the e-expansion naively adding successive orders 
of the perturbation theory series (|J). We show the e 2 — 
and e 3 — approximations. The values predicted by suc- 
cessive orders of the e-expansion for e — 2 do not seem 
to correlate with each other nor with the exact values 
of the d = 2 exponents. It is not only the large expan- 
sion parameter e = 2 that spoils the convergence of the 
perturbation theory. As is well known, the perturbation 
theory series of renormalizcd field theory are asymptotic 
at best Q and have zero radius of convergence |0 . An 
appropriate resummation technique must be used to ex- 
tract reliable data from the series. Here we use a Borel 
resummation refined by a conformal mapping Jl5| as de- 
scribed for this particular case in Ref . . 

The resummation procedure allows us to restore con- 
vergence of the e-expansion and enables us to extract 
reliable exponent values. These are also shown in Fig. [l]. 
The same resummation is applied for the renormalization 
group expansions at fixed D = 2 |Q . All data is summa- 
rized in Table Q. For Af 2 perturbation theory gives the 
exact result at first order with all higher orders vanishing 
Jl8l . For completeness we have included here as in Table 
[n|thc corresponding D = 3 dimensional results from [Q. 

Table || contains the data for the resummed values 
of the exponents A^^ for the combined SAW and RW 
copolymer stars. We show the values obtained by the 
fixed D = 2 technique and the e-expansion in compari- 
son with their exact counterparts (||). Again, the results 
are in fair agreement for not too large numbers of chains 
fa. Indeed as one can see from the formulas (|^)-(^) the 
increase of fa leads to an increase of the e-expansion co- 
efficients resulting in a growing inaccuracy even of the 
resummed series. 

In the perturbative treatment (which starts from the 
upper critical dimension D = 4 j|) the order of the 
chains in the star does not matter for the scaling laws: 
there is always a possibility for every chain to interact 
with any other chain constituting a star. On the other 
hand, in D = 2 the order of chains does matter: each 
chain of the star will interact only with its direct neigh- 
bors. This topological restriction has to be taken into 
account when comparing exact and perturbative results 
for D = 2 linked polymers. 

The data of Tables [jj, [0] convinces us that the perturba- 
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tion theory series for low numbers of chains ji is reliable 
even for D = 2. Also, we note that the D — 2 copoly- 
mers with /i = 2 are worse described by the perturbation 
theory than those with /i = 1 where the topological re- 
striction is not present, due to the symmetry of ordering. 

b. Multifractal Spectrum It is of special interest to 
notice another physical interpretation for the copolymer 
star exponents (Q). For the diffusion of freely moving 
particles in the presence of a polymer absorber it may be 
shown that the moments of particle density display 
universal scaling laws in the vicinity of the absorbing 
polymer. In particular, the nth moment (p n ) at distance 
r from the core of an absorbing polymer star with m 
chains scales as [io|: 



< P » >~ (R/r) 
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here < . . . > denotes the ensemble average over the con- 
figurations of the absorbing polymers while otherwise we 
use the notations of formula ([l]). The harmonic mea- 
sure that corresponds to this absorption phenomenon 
possesses multifractal [|l0[ O, |l2| properties: The flux cf> 
of diffusing particles onto the polymer star is interpreted 
as a measure defined on the fractal structure of the star 
polymer. This flux is proportional to the particle density 
(^) some small distance £ away from an absorbing point 
on the star polymer. The following scaling law for the 
normalized moments of the harmonic measure at the core 
of the star polymer is derived: 



< <j> n > i < <j> > n ~ (R/ey 
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Comparing formulas (||) and (Q) the obvious relation be 
tween the exponents t and A is: 



(8) 



Note, that here we use the standard definition |lO| for the 
spectrum of the exponents r, including normalization of 
the nth moment of the flux in (Q) by < cf> >". The expres- 
sion (0) is a multifractal measure as far as the spectrum r 
is non-trivial (t ^ 0). Obviously, the spectrum A is non- 
trivial as well, as one can see already from the second 
order perturbative results (|j), (||): A mn 7^ nA m i. 

Using the relation (||) we now compare the values of 
the exponents r as obtained in the perturbative approach 
with those of the MC simulations and the exact calcula- 
tions. If the absorbing polymer star is chosen to be RW- 
like, then the spectra r are calculated by substituting the 
exponents A G into (||) (and exponents X u for SAWs, cor- 
respondingly) . We show in Fig. ^| the resulting spectra 
for the spectrum of the exponents for the absorption 
at the end-point of a SAW. Again, we display the exact 
results derived from formulas (J4|) in comparison with the 
perturbative data with and without resummation. We 
observe a similar behavior as for the spectra A: for low 
enough values of n the resummed perturbative renormal- 
ization group results are in reasonable agreement with 
the exact predictions. 



The description of multifractal phenomena often uses 
the spectral function formalism . To obtain this func- 
tion for the absorption process on the center of a star of 
to chains we analytically continue the set of exponents 
T mn in the variable n and calculate the following Legen- 
dre transform 
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According to the standard definition |T(J we have in- 
cluded into ([)]) the fractal dimension d of the absorber. 
In particular, this gives the maximal value of the spec- 
tral function to be equal to d — r m0 = d — A m o- There 
appears to be no natural generalization of d for arbitrary 
to. In our presentation we use that the fractal dimen- 
sion of a polymer star is equal to that of a polymer chain 
(d = 2 for a RW and d = 4/3 for a SAW in D = 2, cor- 
respondingly). Another way to define the multifractal 
spectrum is chosen in Refs. Q where the nth moments 
(Q) are normalized not by < >" but by < <j> >. This 
choice has the remarkable feature that without introduc- 
ing an additional dimension d in (|9|) the maximum of 
f(a) d = 4/3 for both the SAW and the RW absorbers 
with fi — 2. However, this definition is not general, as 
already in D = 3 it does not lead to the correct location 
of the spectral function maximum. 

We display the spectral function for the case of the ab- 
sorption at a RW end in Fig. |j. The spectral function 
is negative for some values of a. This behavior was ob- 
served already in the perturbative studies of Ref. |1C[| 
and related to the ensemble average in formulas (||) , (Q) 
in contrast to the usual site average over the support of 
the measure. Comparing the perturbative, the MC and 
the exact results we note that the region of 1 < n < 5 
in which reasonable agreement for the spectra A and r 
is found corresponds to a rather small region in terms of 
the parameter a in the left wing of the spectral functions. 

c. Conclusions The field theoretical description of a 
polymer network of different species has been derived re- 
cently U in the form of renormalization group pertur- 
bation theory series. In this paper, we compared and ver- 
ified the perturbation expansions for D — 2 with a sub- 
sequently published exact study || and MC simulations. 
We found that the series when appropriately resummed 
give a reliable description of the scaling of copolymer 
stars with not too many chains even for the D = 2 di- 
mensional case. Not only the increase of the expansion 
parameter forbids application of the perturbative results 
for higher numbers of chains: in D = 2 topological re- 
strictions appear that are not taken into account by field 
theory. 

We have analyzed the multifractal properties of the 
resulting spectra of exponents and found that the spec- 
trum r(n) is in much better coincidence with the exact 
counterpart than the spectral function f(a). In the latter 
case the perturbative approach gives comparable results 
only for a narrow region in the left wing of the spectral 
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function. Note however, that the right wing corresponds 
to negative values of n, where obviously the analytic con- 
tinuation of the perturbative expansions as well as of the 
exact results is speculative. 
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h; h 


1;1 1;2 1;3 1;4 1;5 2;2 2;3 2;4 2;5 


exact 


1.25 2 2.693 3.356 4 2.916 3.738 4.510 5.25 


MC 
±A 


1.251 1.986 2.662 3.295 3.908 2.913 3.703 4.506 - 
0.004 0.004 0.005 0.015 0.022 0.005 0.030 0.039 - 


D = 2 

e 3 


1.22 2.00 2.58 3.04 3.43 3.45 4.59 5.52 6.34 
1.20 2.00 2.56 2.99 3.36 3.41 4.49 5.37 6.13 


D = 3 

e 3 


0.58 1.00 1.35 1.69 1.98 1.81 2.53 3.17 3.75 
0.56 1.00 1.33 1.63 1.88 1.77 2.45 3.01 3.51 



TABLE I: Exponents \f lj2 (/i RWs and f 2 RWs) calculated 
by different techniques. Lower part: D = 3 results. 



h; h 


1;1 1;2 1;3 1;4 1;5 2;1 2;2 2;3 2;4 2;5 


exact 


.38 1.04 1.67 2.28 2.88 .88 1.79 2.61 3.39 4.13 


D — 2 

e 3 


0.34 1.01 1.54 2.02 2.41 0.72 1.85 2.80 3.64 4.39 
0.32 1.03 1.57 2.01 2.41 0.76 1.87 2.82 3.62 4.34 


D = 3 

e 3 


0.15 0.51 0.81 1.07 1.32 0.42 1.02 1.57 2.05 2.49 
0.17 0.53 0.81 1.09 1.36 0.42 1.04 1.63 2.15 2.65 



TABLE II: Exponents A£,/ 2 (/i SAWs and / 2 RWs) calcu- 
lated by different techniques. Lower part: D = 3 results. 
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FIG. 1: Exponent A^ n . Lines: exact (bold), bare e 2 
(dashed), and e 3 (dotted) data. Symbols: MC (■, size is 
larger than errorbars); the resummed e 3 (•), and D = 2 data 
(♦) are connected by thin lines. 





FIG. 3: The spectral function for the absorption at a RW 
end. Symbols as in Fig. [[]. 



